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pi I Abstract 

, In this paper we study a subordinate Brownian motion with a Gaussian component and 

I a rather general discontinuous part. The assumption on the subordinator is that its Laplace 

exponent is a complete Bernstein function with a Levy density satisfying a certain growth 
, condition near zero. The main result is a boundary Harnack principle with explicit boundary 

Ph " decay rate for non-negative harmonic functions of the process in C^'^ open sets. As a consequence 

■ of the boundary Harnack principle, we establish sharp two-sided estimates on the Green function 

, of the subordinate Brownian motion in any bounded C ' open set D and identify the Martin 

boundary of D with respect to the subordinate Brownian motion with the Euclidean boundary. 
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1 Introduction 

The infinitesimal generator of a d-dimensional rotationally invariant Levy process is a non-local 



X: 

H. 

. . ■ operator of the form C = 6A -|- A where 6 > and 



Af{x) = / {fix + y)- fix) - Vfix) ■ yl{\y\<i}) vidy) = lim / (/(x + y) - fix)) uidy) . 

The measure v on M'^\{0} is invariant under rotations around origin and satisfies Jj^d (l-^lyP) ^{dy) < 
OO. When = 0, the operator C is proportional to the Laplacian, hence a local operator, while 
when 6 = 0, the operator £ is a purely non-local integro-differential operator. In particular, if 
6 = and vidx) = c\x\~'^~'^dx, a G (0,2), then A is proportional to the fractional Laplacian 
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A"/2 := -(-A)°/2. Levy processes are of intrinsic importance in probability theory, while integro- 
differential operators are important in the theory of partial differential equations. Most of the 
research in the potential theory of Levy processes in the last fifteen years concentrates on purely 
discontinuous Levy processes, such as rotationally invariant stable processes, or equivalently, on 
purely non-local operators of the type A. For summary of some recent results from a probabilistic 
point of view one can consult [U [TOl [Ml ES] and references therein. We refer the readers to [Tj [El [9] 
for a sample of recent progress in the PDE literature, mostly for the case of a fractional Laplacian 
A"/2, a e (0,2). 

In many situations one would like to study operators that have both local and non-local parts. 
From a probabilistic point of view, this corresponds to processes with both a Gaussian component 
and a jump component. The fact that such a process X has both Gaussian and jump components 
is the source of many difficulties in investigating the potential theory of X. The main difficulty 
in studying X stems from the fact that it runs on two different scales: on the small scale the 
diffusion corresponding to the Gaussian part dominates, while on the large scale the jumps take 
over. Another difficulty is encountered when looking at the exit of X from an open set: for 
diffusions, the exit is through the boundary, while for the pure jump processes, typically the exit 
happens by jumping out from the open set. For the process X, both cases will occur which makes 
the process X much more difficult to study. 

Despite the difficulties mentioned above, in the last few years significant progress has been made 
in understanding the potential theory of such processes. Green function estimates (for the whole 
space) and the Harnack inequality for a class of processes with both diffusion and jump components 
were established in [27\ I31| . The parabolic Harnack inequality and heat kernel estimates were 
studied in [32] for Levy processes in M'^ that are independent sums of Brownian motions and 
symmetric stable processes, and in [16] for much more general symmetric diffusions with jumps. 
Moreover, an a priori Holder estimate was established in [16] for bounded parabolic functions. For 
earlier results on second order integro-differential operators, one can see [17] and the references 
therein. 

Important progress has been made in two recent papers [121 [13] which consider operators of the 
type A -|- a"A"/^ for a G [0, M]. The process corresponding to such an operator is an independent 
sum of a Brownian motion and a rotationally invariant a-stable process with weight a. In [12] the 
authors established a (uniform in a) boundary Harnack principle (BHP) with explicit boundary 
decay rate for non-negative harmonic functions with respect to A -|- a°A°/2 in C ' open sets. 
By using the BHP, the second paper [13] established sharp Green function estimates in bounded 
C^'^ open sets D, and identified the Martin boundary of D for the operator A -|- a"A"/^ with its 
Euclidean boundary. 

The purpose of the current paper is to extend the results in [El [13] to more general operators 
than A-|-a"A"/^. Analytically, the operators that we consider are certain functions of the Laplacian. 
To be more precise, we consider a Bernstein function (p : (0, oo) — )• (0, oo) with 4>{0+) = 0, i.e., is 
of the form 

^{X) = b\+ I {1 - e-^^) fi{dt) , A>0, (1.1) 

i(0,oo) 

where 6 > and /i is a measure on (0, oo) satisfying /(ooo)(^ ^ f^i^t) < oo. /i is called the Levy 
measure of (p. By Bochner's functional calculus one can define the operator (p{A) := —(/>(— A) which 
on C^iJR.'^), the collection of bounded functions in with bounded derivatives, turns out to be 
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an integro-differential operator of the type 

bAfix) + [ {fix + y)- fix) - Vfix) ■ 2/l{|,|<i}) ^dy) , 
where the measure ly has the form vidy) = ji\y\) dy with j : (0, oo) — t- (0, oo) given by 

/>oo 

jir) = / (47rt)-^/2g-rV(4t) _ 
Jo 

In order for the operator to have both local and non-local parts we will assume that 6 > and 
fi ^ 0. In fact, without loss of generality, throughout the paper we always suppose that 6=1. 
Note that by taking (/)(A) = A + a"A"/^ we are back to the operator A + a"A"/^. 

The operator i;^(A) is the infinitesimal generator of the Levy process X that can be constructed 
as follows: Recall that a one-dimensional Levy process S = iSt : t > 0) is called a subordinator if 
it is non-negative and = 0. A subordinator 5 can be characterized by its Laplace exponent (p 
through the equality 

E[e-^^'] = e-*'^(^), t>0,X>0. 

The Laplace exponent (p can be written in the form (jl.ip . We will assume that 6=1. Suppose 
that W = iWt : i > 0) is a d-dimensional Brownian motion and S = (S'j : i > 0) is a subordinator, 
independent of W, with Laplace exponent (j). The process X = {Xt : t > 0) defined by Xt = Ws^ 
is called a subordinate Brownian motion and its infinitesimal generator is </>(A). It is a sum of a 
Brownian motion and an independent purely discontinuous (rotationally invariant) Levy process. 

Potential theory of one-dimensional subordinate Brownian motions in this setting was studied 
in [23]. In the current paper we look at the case when d > 2. In order for our methods to work 
we need additional assumptions on the Bernstein function (p. We will assume that is a complete 
Bernstein function, namely that the Levy measure fi has a completely monotone density. By a 
slight abuse of notation we will denote the density by For the Levy density // we assume a 

growth condition near zero: For any K > 0, there exists c = c{K) > 1 such that 

/i(r) < c^i(2r), VrG(0,Er). (1.2) 

We will later explain the role of these additional assumptions. 

To state our main result, we first recall that an open set D in (when d > 2) is said to be C^'^ 
if there exist a localization radius i? > and a constant A > such that for every Q G dD, there 
exist a C^'^-function ip = ipg : R'^-i R satisfying ip{0) = 0, V(^(0) = (0, ...0), \\Vip\\oo < A, 
\V^p{x) — Vip{y)\ < A\x — y\, and an orthonormal coordinate system CSq: y = {yi, ■ ■ ■ ,yd-i,yd) ='■ 
(y^ yd) with its origin at Q such that 

B{Q, R)nD = {y= (y, y^) G S(0, R) in CSq : y^ > ip{y)}- 

The pair (R, A) is called the characteristics of the C^'^ open set D. Note that a C^'^ open set can 
be unbounded and disconnected. 

For any x € D, 5d{x) denotes the Euclidean distance between x and W^. For any x ^ D, 
5q£){x) denotes the Euclidean distance between x and dD. It is well known that, if Z) is a C^'^ 
open set D with characteristics {R,A), there exists R< R such that D satisfies both the uniform 
interior ball condition and the uniform exterior ball condition with radius R: for every x £ D with 
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^d{x) < R and y G M°' \ Z) with 6QD{y) < R, there are Zx,Zy S dD so that \x — Zx\ = Sd{x), 
\y — Zy\ = SgD{y) and that B{xo, R) C D and B{yo,R) C W^\D where xq = Zx + R{x — Zx)/\x — Zx\ 
and yo = Zy + R{y — Zy)/\y — Zy\. Without loss of generahty, throughout this paper, we assume 
that the characteristics (i?, A) of a C^'^ open set satisfies R = R < 1 and A > 1. 

For any open set D C M.'^, td := inf{t > : Xt ^ D} denotes the first exit time from D by X. 

Definition 1.1 A function / : R'^ i— t- [0, oo) is said to be 

(1) harmonic in an open set D C M"^ with respect to X if for every open set B whose closure is a 

compact subset of D, 

/(x) =E,[/(X,^)] for every xeB- (1.3) 

(2) regular harmonic in D for X if for each x £ D, f{x) = [/(X^-^)]. 

We note that, by the strong Markov property of X, every regular harmonic function is automatically 
harmonic. 

Let Q G dD. We will say that a function / : M'^ — )■ M vanishes continuously on D'^ n B{Q, r) if 
/ = on n B{Q,r) and / is continuous at every point of dD n B{Q,r). The following is the 
main result of this paper. 

Theorem 1.2 Suppose that the Laplace exponent (j) of the subordinator S, independent of the 
Brownian motion W , is a complete Bernstein function and that the Levy density of S satisfies 
(jl.2p . Let X = {Xf : t > 0) be the subordinate Brownian motion defined by Xt = W{St). For any 
C^'^ open set D in M.'^ with characteristics {R,A), there exists a positive constant C = C{d,A,R) 
such that for r G {0,R\, Q G dD and any nonnegative function f in which is harmonic in 
D n B{Q,r) with respect to X and vanishes continuously on D'^ n B{Q,r), we have 

fix) ^(.py) for every x,y e DnB{Q,r/2). (1.4) 



Snix) Sniy) 

We note that (jl.4p is a strengthened version of the usual boundary Harnack principle stated for 
the ratio of two non-negative functions, / and g, harmonic in D n B{Q,r) with respect to X, and 
which says that 

f(x) f(v) 

jy)^(jJ_^ for every x,y £ DnB(Q,r/2). 
9{x) g{y) 

Indeed, the above inequality is a consequence of (II. 4|) . We note that (jl.4|) gives the precise boundary 
decay of non-negative harmonic functions and that the function x i— >• Soix) is not harmonic in 
D n B{Q, r) with respect to X. 

Remark 1.3 The same type of boundary Harnack principle in C^'^ domains is valid also for 
Brownian motions, namely the boundary decay rate is of the order 6d{x)- Since on the small scale 
the diffusion part of X dominates, one would expect that harmonic functions of X and of Brownian 
motion have the same decay rate at the boundary. For this reason, some people might expect that 
some kind of perturbation methods can be used to prove the BHP for X. We note that it is unlikely 
that any perturbation method would work because of the following: Suppose that instead of X we 
consider a process X" with the infinitesimal generator 

Cy{x) = A/(x) + / {fix + y)- fix) - Vfix) ■ yl{\y\<i}) v'^idy) , 
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where i'"'{dy) = l||y|<a} with < a < oo. Thus X"' is the process X with jumps of size 

larger than a suppressed. In Section [6] we present an example of a (bounded) C^'^ domain D on 
which the boundary Harnack principle for X°' fails, even for regular harmonic functions vanishing 
on D^. Note that if we think of X as a perturbation of Brownian motion, then X"" is an even 
smaller perturbation of the same Brownian motion. The counterexample in Section [6] shows that, 
despite the (seemingly) local nature of the BHP, one needs some information of the structure of 
large jumps of X. 

For any open set D C M*^, we will use X^ to denote the process defined by (w) = Xt{uo) if 
t < T£,[uj) and XP{uj) = d if t > T£){uj), where 9 is a cemetery point. The Green function of X^ 
will be denoted by G]^{x,y). For the precise definition of Gd, see Section 2. 

To state our result on Green function estimates, we introduce a function go first. For d > 2, 
we define for x,y G D, 

f ( 1 A ^af^OMi) \ ^hen d>3, 

an(x v) = I 1^-4 V l^-^l / 
' ^ log l + whend = 2. 



Theorem 1.4 Suppose that the Laplace exponent (p of S is a complete Bernstein function and 
that the Levy density of S satisfies (|1.2p . For any bounded C^'^ open set D C W^, there exists 
C = C{D) > 1 such that for all x,y £ D 

C~^gD{x,y) < GD{x,y) < C goix^y). (1.5) 

Finally, we state the result about the Martin boundary of a bounded C^'^ open set D with 
respect to X^ . Fix xq £ D and define 

Moix.y) := x,yeD,yj^xo. 

A function / is called a harmonic function for X^ if it is harmonic for X in D and vanishes outside 
D. A positive harmonic function / for X^ is minimal if, whenever g is a positive harmonic function 
for X^ with g < f on D, one must have f = eg for some constant c. 

Theorem 1.5 Suppose that D is a bounded C^'^ open set in W^. For every z € dD, there exists 
Md{x, z) := Ymiy^z M£){x, y). Further, for every z € dD, M/)(-, z) is a minimal harmonic function 
for X^ and M£){-,zi) 7^ M/)(-,zi) if zi 7^ Z2- Thus the minimal Martin boundary of D can be 
identified with the Euclidean boundary. 

Thus, by the general theory of Martin representation in j26j and Theorem 1 1.5 1 we conclude that, 
for every harmonic function n > with respect to X^, there is a unique finite measure u on dD 
such that u{x) = Jq^^ M£){x, z)i'{dz). 

Let us now describe the main ingredients of the proof of Theorem 11.21 the boundary Harnack 
principle. We follow the general strategy for proving the boundary Harnack principle in different 
settings which requires the Carleson estimate, and upper and lower estimates on exit probabilities 
and exit times from certain sets usually called "boxes" (see \5\ [T2 \ \18 \ I20|) . In Theorem 15.31 we 



5 



prove the Carleson estimate for a Lipschitz open set by modifying the proof in [12]. In order to 
obtain the upper exit probabihty and exit times estimates, we fohow the approach from |12] . the 
so-called "test function" method (which was modeled after some earlier ideas, see [UlIS]), but have 
to make major modifications. In [12], the test functions are power functions of the form x i— )■ {x^)^ 
which are either superharmonic or subharmonic for the corresponding process, and the values of 
the generator on these test functions are computed in detail. In our setting, the power functions are 
neither superharmonic nor subharmonic, and explicit calculations cannot be carried out because of 
the lack of explicit form of the Levy measure. Instead we use the approach developed in |25j for 
the case of certain pure-jump subordinate Brownian motions, which seems to be quite versatile to 
cover various other cases. 

One of the main ingredients in [25j comes from the fluctuation theory of one-dimensional Levy 
processes. Its purpose is to identify a correct boundary decay rate by finding an appropriate 
harmonic function. Let Z = {Zt ■ t > 0) be the one-dimensional subordinate Brownian motion 
defined by Zt := W^^, and let V be the renewal function of the ladder height process of Z. The 
function V is harmonic for the process Z killed upon exiting (0,oo), and the function w{x) := 
V{xci)l^r^^^Q'^, X G M'^, is harmonic for the process X killed upon exiting the half space M^J. := {x = 
{xi, . . . ,Xci-i,Xd) G M*^ : Xrf > 0} (Theorem 13. 2p . Therefore, w gives the correct rate of decay of 
harmonic functions near the boundary of M^J.. We will use the function w as our test function. 
Note that the assumption that cp is a complete Bernstein function implies that w is smooth. Using 
smoothness and harmonicity of w together with the characterization of harmonic functions recently 
established in |llj . we show that (A + A)w = on the half space (Theorem 13. 4p . Consequently we 
prove the following fact in Lemma |4.H which is the key to proving upper estimates: If D is a C^'^ 
open set with characteristics {R,A), Q G dD and h{y) = V{5Diy))iDnB{Q,R)j then (A + A)h{y) is 
a.e. well defined and bounded for y £ D close enough to the boundary point Q. Using this lemma, 
we give necessary exit distribution estimates in Lemma 14.31 Here we modify the test function h 
by adding a quadratic type function (in one variable) - this is necessary due to the presence of 
the Laplacian. The desired exit distribution estimates are directly derived by applying Dynkin's 
formula to the new test function. The reader will note that our proof is even shorter than the one 
in [12], partly because, in |12j . the uniformity of the boundary Harnack principle for A + a"A"/^ 
in the weight a G (0, M] is established. 

In order to prove the lower bound for the exit probabilities we compare the process X killed upon 
exiting a certain box D with the so-called subordinate killed Brownian motion obtained by first 
killing Brownian motion upon exiting the box D, and then by subordinating the obtained process. 
If the latter process is denoted by , then its infinitesimal generator is equal to — (?!)(— A|^). Here 
A|^ is the Dirichlet Laplacian and —(j){—A\^) is constructed by Bochner's subordination. The 
advantage of this approach is that the exit probabilities of X^ dominate from the above those of 
the process Y^, while the latter can be rather easily computed, see [33]. This idea is carried out 
in Lemma 14.41 (as well as for some other lower bounds throughout the paper). 

Once the boundary Harnack principle has been established, proofs of Theorems 11.41 and 1 1 . 5 1 are 
similar to the corresponding proofs in |13] for the operator A + a" A". Therefore we do not give 
complete proofs of these two theorems in this paper, only indicate the necessary changes to the 
proofs in [13]. 

The rest of the paper is organized as follows. In the next section we precisely describe the 
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settings and recall necessary preliminary results. Section 3 is devoted to the analysis of the process 
and harmonic functions in the half-space. Section 4 is on the analysis in C^'^ open sets, and is 
central to the paper, and this is where most of the new ideas appear. In this rather technical section 
we establish the upper and lower bounds on the exit probabilities and exit times. In Section 5 we 
first prove the Carleson estimate for Lipschitz open sets and then prove the main Theorem 11.21 In 
Section 6 we provide the counterexample already mentioned in Remark [6l Finally, in Section 7 we 
explain the differences between the proofs of Theorems 11.41 and 11.51 and their counterparts from 

Throughout this paper, the constants Ci, C2, R, Ri, R2, R3 will be fixed. The lowercase 
constants ci,C2,--- will denote generic constants whose exact values are not important and can 
change from one appearance to another. The dependence of the lower case constants on the 
dimension d may not be mentioned explicitly. We will use ":=" to denote a definition, which is 
read as "is defined to be". For a, 5 € M, a A 6 := min{a,6} and a V 6 := max{a, 5}. For every 
function /, let /"*" := / V 0. For every function /, we extend its definition to the cemetery point 
d by setting f{d) = 0. We will use dx to denote the Lebesgue measure in and, for a Borel set 
A C M.'^, we also use \A\ to denote its Lebesgue measure. 

2 Setting and Preliminary Results 

A C°° function (p : (0, 00) — t- [0,oo) is called a Bernstein function if (— 1)"L'"'<^ < for every 
n = 1, 2, ... . Every Bernstein function has a representation 4>{X) = a + bX + Jj-q ^^-((l — e""^*) fi{dt) 
where a, 6 > and ^ is a measure on (0, 00) satisfying Jj-q ^^^(l At) fi{dt) < 00; a is called the killing 
coefficient, b the drift and fi the Levy measure of the Bernstein function. A Bernstein function (p is 
called a complete Bernstein function if the Levy measure fi has a completely monotone density 
i.e., (— l)"'L'"/i(t) > for every non- negative integer n and all t > 0. Here and below, by abuse 
of notation we denote the Levy density by fJ.{t). For more on Bernstein and complete Bernstein 
functions we refer the readers to |29j . 

A Bernstein function (p on (0, 00) is the Laplace exponent of a subordinator if and only if 
4>{0+) = 0. Suppose that S is a subordinator with Laplace exponent (p. S is called a complete 
subordinator if is a complete Bernstein function. The potential measure U of S is defined by 

POD 

U{A)=K lis,eA}dt, ^C[0,oo). (2.1) 
Jo 

Note that U{A) is the expected time the subordinator S spends in the set A. 

Throughout the remainder of this paper, we assume that S = {St : t > 0) is a complete 
subordinator with a positive drift and, without loss of generality, we shall assume that the drift of 
S is equal to 1. Thus the Laplace exponent of S can be written as 

0(A):=A + ^(A) where V'(A) := / (1 - e'^^) n{dt). 

i{0,oo) 

We will exclude the trivial case of St = t, that is the case of "0 = 0. Since the drift of S is equal 
to 1, the potential measure U of S has a completely monotone density u (cf. [6l Corollary 5.4 and 
Corollary 5.5]). 
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Suppose that W = {Wt : t > 0) is a Brownian motion in M'^ independent of S and with 

E,.[e^^-(^'-'^o)] = e-*l^l', for all x,e e M'^. 

The process X = (Xt t > 0) defined by Xt = Wst is called a subordinate Brownian motion. It 
follows from [6l Chapter 5] that X is a Levy process with Levy exponent = |0p + "0(1^1^): 

^^l^ie-ix^-Xo)^ = e-^-^d^l'), for aU x, 9 G R". 

The Levy measure of the process X has a density J, called the Levy density, given by J{x) = j{\x\) 
where 



poo 

:= / (47rt)- 



■d/2-ry{4t) 



fi{t) dt, 



r > 0. 



(2.2) 



Note that the function r i— )• j{r) is continuous and decreasing on (0,oo). We will sometimes use 
the notation J{x,y) for J{x — y). 

The function J(x, y) is the Levy intensity of X. It determines a Levy system for X, which 
describes the jumps of the process X: For any non-negative measurable function / on M+ X M*^ X M"' 
with /(s, y,y) =0 for all y G M'^, any stopping time T (with respect to the filtration of X) and any 
X G M'^, 

cT 



s<T 



E, 







fis,Xs,y)J{Xs,y)dy ds 



(2.3) 



(See, for example, [21 Proof of Lemma 4.7] and [151 Appendix A].) 

Recall that for any open set U C W^, tjj = inf{t > : Xt ^ is the first exit time from U by 
X. The following simple result will be used in Section 5. 

Lemma 2.1 For every q > Q, there exists c = c{g) > such that for every xq G M'^ and r G (0, g], 



1^2 



c r 



< [TB{xo,r)] < C) 



(2.4) 



Proof. In the case d> 3, this lemma has been proved in ^7\. Moreover, one can easily adapt the 
proofs of j3H Lemmas 2.1-2.1] to arrive at the desired lower bound for all dimensions. Here we 
provide a proof of the desired upper bound that works for all dimensions. 

Let Cq{W^) be the collection of functions in M*^ vanishing at infinity. For any / G Cq(M'^), 
we define 

Af{y)= I' {f{z + y)-f{y)-z-Vf{y)l{\,\^i})J{z)dz, yeR''. 
Let > 5 be such that 

/ \z\'^J{z)dz < -N^. 
J{\z\<i} 2 

Let 5 be a radial function taking values in [0, 2] such that 




\y\ < 1 
2 < \y\ < 3 
\y\> N - 1. 
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For any r > 0, put f{y) = g{y/r). Then for y G 5(0, r), Af{y) = 2dr ^. For any y S i?(0, r), we 
have 



\Af{y)\ = \ {f{z + y)- f{y)-z-Vf{y)l{\,\^Nr})J{z)dz\ 



< \ ^f(^z + y)-fiy)-z-Vfiy))Jiz)dz\+ / f{y)J{z)dz 

J{\z\<Nr} J{Nr<\z\} 

< cxr""^ / |zpJ(z)dz + r-2|y|2 / J[z)dz 

J{\z\<Nt} J{Nr<\zW 

< cir-2 [ \z\'^J{z)dz + N-^r^^ [ \z\'^J{z)dz+ [ J{z)dz. 

J{\z\<Nr} J{Nr<\z\<l} •^{|z|>l} 

Thus we know that there exist ro G (0, 1) and C2 > such that for any r € (0, tq), 

A/(y) + ^/(y) > C2r-2, yGBiO,r). 

Using this and the fact that A + ^ is the infinitesimal generator of the process X, by the Dynkin's 
formula, we have that for r £ (0, tq), 

^xo['rB{xo,r)] = IEo[rB(o,r)] = limlEo[rij(o,r) At] 

< c^V^limEo / {A + A)f{Xs)ds 

ttoo Jq 

tJOQ 

Now the desired upper bound follows easily. □ 

In the remainder of this paper, we will need some control on the behavior of j near the origin. 
For this, we will assume that for any K > 0, there exists c = c{K) > 1 such that 

H{r) <cfj,{2r), VrG(0,i^). (2.5) 

On the other hand, since i;^> is a complete Bernstein function, it follows from \25\ Lemma 2.1] that 
there exists c > 1 such that /i(t) < c/i(t + 1) for every t > 1. Thus by repeating the proof of 
[271 Lemma 4.2] (see also [Ml Proposition 1.3.5]), we can show that for any K > 0, there exists 
c = c{K) > 1 such that 

j{r)<cj{2r), VrG(0,K), (2.6) 

and, there exists c > 1 such that 

j{r) <cj{r + l), Vr>l. (2.7) 

Note that, as a consequence of (j2.6p . we have that, for any K > 0, 

j{ar) < ca'^jir), Vr G {0,K) and a G (0, 1) (2.8) 

where c = c{K) is the constant in (j2.6p and u = ^{K) := log2 c. 

The following Harnack inequality will be used to prove the main result of this paper. 
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Proposition 2.2 (Harnack inequality) There exists a constant c > such that for anyr G (0, 1] 
and xq G M"^ and any function f which is nonnegative in R"' and harmonic in B{xQ,r) with respect 
to X we have 

f{x) < cf{y) for all x,y e B{xq, r/2). 

Proof. We first deal with the case (i > 3. When / is bounded, this proposition is just |271 Theorem 
4.5]. Using the same argument as in the proof of \n\ Corollary 4.7], one can easily see that \n\ 
Theorem 4.5] can be extended to any nonnegative harmonic function. 

The assertions of the proposition in the cases of d = 2 and d = 1 follow easily from the 
assertion in the case d > 3. Since the arguments are similar, we will only spell out the details in 
the case d = 2. For any x G M^, x = {x,x'^), where x G M^. Analog notation will be used also 
for other objects in M^. Let X = {Xt,'Px) be the subordinate Brownian motion in and write 
X = {X,X^). Note that X has the same distribution under IP(s^o) ^'^d '^(x^x^) for any G M. 
Hence we can define := P(x,o)- The process (X,Pj) is a subordinate Brownian motion in via 
the same subordinator as the one used to define X. For any given / : — t- [0, oo), we extend it to 

by defining fix) = f{{x,x^)) := fix). Then 

(1) If / is regular harmonic (with respect to X) in an open set D C M^, then / is regular harmonic 

(with respect to X) in the cylinder D := D >iM.. Indeed, let := inf{i > : Xt ^ D} he 
the exit time of X from D, and td := inf{t > : Xt ^ D}. Then clearly = T£,. Thus, for 
any x = (x, x^) G D, 

^x[f(Xr,)] = ^x[f{Xr~)] = m = fix) . 

(2) If / is harmonic (with respect to X) in an open set D cM?, then / is harmonic (with respect 

to X) in the cylinder D := D x R. Indeed, let i? C .D be open and relatively compact. Then 
there exists a cylinder C = C x M such that B <Z C and C C .D is open and relatively compact 
(in D). Since / is harmonic (with respect to X) in D, it is regular harmonic in C. By (1), / 
is regular harmonic (with respect to X) in C, and therefore also harmonic in C. Since B is 
compactly contained in C, we see that 

/(x) =E,[/(X,^)], foraUxGS. 

Let r G (0,1), xq £ and define xq := (xO)0). Assume that / : — )• [0, oo) is harmonic 
(with respect to X) in BixQ,r). Then / defined by /(x) = /(x) is harmonic in B{xQ,r) x R. In 
particular, / is harmonic in B{xQ,r). By the assertion in the case d = 3, 

f{x) < cf{y) , for all x, y G B{xo, r/2) . 

Let X, y G B{xo,r /2), and define x := (x, 0), y := (y, 0). Then 

fix) = fix) < cfiy) = cfiy) . 

□ 

It follows from [U Chapter 5] that the process X has a transition density pit,x,y), which is 
jointly continuous. Using this and the strong Markov property, one can easily check that 

PDit,x,y) :=pit,x,y) - Ex[pit - td, Xrjy,y);t > td], x,y G D 
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is continuous and is the transition density of . For any bounded open set D C M"', we will use 
to denote the Green function of , i.e., 



Note that Gd{x, y) is continuous on {(x, y) £ D x D : x ^ y}. 

3 Analysis on half-space 

Recall that X = (Xt : t > 0) is the (i-dimensional subordinate Brownian motion defined by 
Xt = Wst, where W = {W^, . . . , W^) is a d-dimensional Brownian motion and S = {St : i > 0) an 
independent complete subordinator whose drift is equal to 1 and whose Levy density satisfies ()1.2p . 

Let Z = {Zt : t > 0) be the one-dimensional subordinate Brownian motion defined as Zt := W^^. 
Let Zt := sup{0 V : < s < t} be the supremum process of Z and let L = (Lt : i > 0) be a 
local time of Z — Z at 0. L is also called a local time of the process Z reflected at the supremum. 
The right continuous inverse L^^ of L is a subordinator and is called the ladder time process of 
Z. The process Ht = Z ^-i is also a subordinator and is called the ladder height process of Z. 
(For the basic properties of the ladder time and ladder height processes, we refer our readers to [H 
Chapter 6].) The ladder height process H has a drift i \2.2>\ Lemma 2.1]). The potential measure of 
the subordinator H will be denoted by V. Let V{t) := V{{0,t)) be the renewal function of H. 

By [U Theorem 5, page 79] and |23[ Lemma 2.1], V is absolutely continuous and has a continuous 
and strictly positive density v such that v{0+) = 1. The functions V and v enjoy the following 
estimates near the origin. 

Lemma 3.1 ( |23| Lemma 2.2]) Let R > 0. There exists a constant c = c{R) > 1 such that for all 
X £ (0, R], we have < v{x) < c and c~^x < V{x) < cx . 

By [251 Proposition 2.4], the Laplace exponent x of the ladder height process H of Zt is also a 
complete Bernstein function. Using this and the fact that x has a drift, we see from [241 Corollary 
2.3], that V is completely monotone. In particular, v and the renewal function V are C°° functions. 

We will use to denote the half-space {x = (xi, . . . , x^-i, x^) := {x,Xd) E M'^ : > 0}. 
Define w{x) := y((x^)+). 

Theorem 3.2 The function w is harmonic in with respect to X and, for any r > 0, regular 
harmonic in W^~^ x (0,r) with respect to X. 

Proof. Since Zt := has a transition density, it satisfies the condition ACC in [30], namely 
the resolvent kernels are absolutely continuous. The assumption in [30] that is regular for (0, oo) 
is also satisfied since X is of unbounded variation. Further, by symmetry of Z, the notions of 
coharmonic and harmonic functions coincide. Now the theorem follows by the same argument as 
in [251 Theorem 4.1]. O 

Unlike [251 Proposition 4.2], we prove the next result without using the boundary Harnack 
principle. 




x,y £ D. 
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Proposition 3.3 For all positive constants tq and L, we have 



sup 

0<Xa<L J B{x,roy 



/ w{y)j{\x -y\)dy 

J Blx.mYrX.'^, 



< oo . 



Proof. Without loss of generality, we assume x = 0. We consider two separate cases, 
(a) Suppose L > > ro/4. By ()2.3p and Theorem 13.21 for every x G M^|_, 



w{x) > 
= 



w(Xr 



B(x,rQ/2)nM.\ 



) -.Xr , G B{x,roy 







j{\Xt - y\)w{y) dydt 



(3.1) 



Since 1^ — y| < |x — z| + |x — i/| < + \x — y\ < 2\x — y\ for {z, y) G B{x, ro/2) x B{x, vqY, using 
(j2.6p and (j2.7p . we have j{\z — y\) > cij{\x — y\). Thus, combining this with (j3.ip . we obtain that 



B{x,ro)''nR'l 



'w{y)j{\x - y\)dy < q ^ 



w[x) 



V{L) 



['^ B{x,ro/2)rm.l ] B(0,ro/4)] 



(b) Suppose Xd < ro/4. Note that if |y — x| > tq, then \y\ > \y — x\ — \x\ > 3ro/4 and \y\ < \y — x\ + 
\x\ < |y — x|+ro/4< \y — x\ + \y — x\/4:. Thus, using (j2.6p and (j2.7p . we have j(|y|) > C2j(|3^ ~ vl) 
and 



sup / w{y)j{\x - y\)dy < 03 w{y)j{\y\)dy. 

0<Xd<ro/4 JB(i:,ro)'=nM!{. JB(0,ro/2)=nIR^ 



(3.2) 



Let xi := (0,ro/8). By Theorem O and (fTS]) . 



oo > tt;(xi) > E^.^ 
= E, 



,):Xr , e B(x,ro/2f n: 



'''B(0,ro/4)nMl 



S(0,ro/2)'=n]Rj 



j'd^t - y\)w{y)dydt 



(3.3) 



Since |z - y| < l^] + |y| < (ro/4) + \y\ < 2\y\ for {z,y) G 5(0, ro/4) x 5(0,ro/2)^ using and 
(j2.7p . we have j(k ~ y|) ^ csjdyl). Thus, combining this with ()3.3p . we obtain that 



oo > t(;(xi) > caE^.^ 

= C3E^Jr^(o_^^ 
Combining (j3.2p and (j3.4p . we conclude that 



S(0,ro/4)nM° 



B(0,ro/2)=nIRj 
B(0,ro/2)=nlR^ 



sup /" ?i'(y)i(k - y\)dy < ^(^o/^) ^ < 



(3.4) 



a;eM<*: 0<a;d<ro/4 JB{x,roy 



IEo[TB(0,ro/8)] 



□ 



12 



We now define an operator (A + 2)(A + A)) as follows: 
Af{x) := lim / (/(y) - /(x)) j{\y - x\) dy, 

JBix.eY 



?(x,e)<= 

S)(A + ^) := J / G C2(M^) : lim 



/ (/(y) ~ /(^)) ~ A) '^y exists and is finite > . (3.5) 

JB{x,eY J 



Recall that Cq{W^) is the collection of functions in vanishing at infinity It is well known 
that Cl(W'-) C D(A + ^) and that, by the rotational symmetry of X, A + ^ restricted to Cl{W^) 
coincides with the infinitesimal generator of the process X (e.g. [281 Theorem 31.5]). 

The proof of the next result is similar to that of [251 Theorem 4.3]. We give the proof here for 
completeness. 

Theorem 3.4 (A + A)w{x) is well defined and (A + A)w{x) = for all x G M^. 

Proof. It follows from Proposition 13.31 and the fact that j is a Levy density that for any L > 
and e G (0,1/2) 



sup 



{w{y) -w{x))j{\y - x\)dy 

B{x,eY 



< sup / w{y)j{\y - x\)dy + V{L) j{\y\)dy <oo. (3.6) 

xPR'^: 0<x^<L-JB(x,eY J B(x,eY 



0<Xi<LJ B{x,eY J B{x,e) 

Hence, for every e G (0,1/2), A£w{x) := Jb{x sY^^^^^ ~ ~ x\)dy is well defined. Using 

the smoothness of w in and following the same argument in \25\ Theorem 4.3], we can show 
that Aw is well defined in and A£'w{x) converges to 

Aw{x) = / {w{y) - w{x) - l{|j^„^|<i}(y - x) ■ Vw{x)) j{\y - x\)dy 

locally uniformly in W!^ as e — )• and the function Aw{x) is continuous in M!^. 

Suppose that Ui and U2 are relatively compact open subsets of MJ^ such that Ui C U2 C U2 C 
M^i_. It follows again from the same argument in [25^ Theorem 4.3] that the conditions [111 (2.4), 
(2.6)] are true. Thus, by [111 Lemma 2.3, Theorem 2.11(ii)], we have that for any / G C^(]R!^), 

0=/ Vw{x)-Vf{x)dx + l [ [ iw{y)-wix)){f{y)- f{x))j{\y-x\)dxdy. (3.7) 

For / G C2(]R^) with supp(/) C C f/2 C C M^, 

\w{y) - w{x)\\f(%)) - /(x)lj(l2/ - x\)dxdy 

\w{y) - 7i;(x)ll/(y) - /(x)lj(l?/ - x\)dxdy + 2 / \w{y) - w{x)\\f{x)\j{\y - x\)dxdy 
U2JU2 JuiJui 

<_ej |,-.5-(|,-.|)...,.2||/|U|r,,|(™p„„.))/„|,-.|)., 

JU2XU2 \xeUi J Ju^ 

+ 2II/II00 / / w{y)3{\x - y\)dydx 
Jui Ju^ 
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is finite by Proposition 13.31 and the fact that j{\x\)dx is a Levy measure. Thus by (j3.7p . Fubini's 



theorem and the dominated convergence theorem, we have for any / E c}'^'^ 



0=1 Vw{x) ■Vf{x)dx + l-lim [ {w{y) - w{x)){f{y) - f{x))j{\y - x\) dx dy 

JR'' ^ ^-1-0 J{(x,y)m'^xR'i, \y-x\>£} 



Aw{x)f{x) dx — hm / f{x) / {w{y) — w{x))j{\y — x\)dy dx 

^■^^ jRj \JB{x,eY J 

/ Aw{x)f{x)dx— / f{x)Aw{x)dx = — {A + A)w{x)f{x) dx 

jRd J^d J^d 



where we have used the fact A^w — >■ Aw converges uniformly on the support of /. Hence, by the 

id 



continuity of (A + A)w, we have (A + A)w{x) = in R'^. □ 



4 Analysis on C open set 

Recall that A > 1 and that D is a C^'^ open set with characteristics {R, A) and D satisfies the 
uniform interior ball condition and the uniform exterior ball condition with radius R < 1. The 
proof of the next lemma is motivated by that of [25\ Lemma 4.4]. 

Lemma 4.1 Fix Q G dD and define 

h{y) := V (doiy)) lDni?(0,/?)(y)- 

There exists Ci = Ci{A,R) > independent of Q such that (A + A)h is well defined in D D 
B{Q,R/A) a.e. and 

\{A + A)h{x)\ < Ci for a.e. x e DnB{Q,R/A). (4.1) 

Proof. In this proof, we fix x € D Pi B{Q, -R/4) and xq £ dD satisfying 6d{x) = \x — xo\. We also 
fix the C^'^ function if and the coordinate system CS = CSxq in the definition of C^'^ open set so 
that x = {0,Xd) with < xa < R/4: and B{xo,R)nD = {y = (y, yd) G B{0,R) in CS : yd> V^{y)}- 
Let 

(^i(y) :=i?- V^2_|y|2 and ip2{y) := -R + R^ - \y\^ . 

Due to the uniform interior ball condition and the uniform exterior ball condition with radius R, 
we have 

My) < ^(y) < My) for every y£Dn B{x,R/A). (4.2) 
Define H+ := {y = {y, yd) £CS:yd>0} and let 

A:={y = (y, yd) £ (D U H+) n B{x, R/A) : My) < yd < Mv)}, 

E:={y = {y,yd) G B{x,R/A) : yd > My)}- 

Note that, since \y - Q\ < \y - x\ + \x - Q\ < R/2 for y G B{x, R/4:), we have B{x, R/A) D D C 
B{Q,R/2)r\D. 
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Let 

h.{y) ■.= V[5^^{y)). 

Note that hx{x) = h{x). Moreover, since 6^^{y) = {yd)~^ in CS, it follows from Theorem 13.41 that 
Ahx is well defined in and 

(A + ^)/i,(y) = 0, yy&H+. (4.3) 
We show now that A{h — hx){x) is well defined. For each e > we have that 



{y(iD\jH+:\y-x\>e} 



{h{y) - hx{y))j{\y - x\) dy 



< / (Hy) + hxiy))j{\y - x\)dy + / {h{y) + hxiy))j{\y - x\) dy 

JB{x,R/AY J A 

+ / \h{y) - hx{y)\j{\y - x\)dy =: h + h + h- 
Je 

By the fact that h{y) = for y € B{Q, Rf, 

h< sup / V{yd)j{\z-y\)dy + ci j{\y\)dy =: Ki + K2. 

K2 is clearly finite since J is the Levy density of X while Ki is finite by Proposition 
For y £ A, since V is increasing and {R — \/ R"^ — < R^^\y\'^, we see that 

K{y) + h{y) < 2V[My) - My)) < 2Vi2R-^\y\^) < 2V{2R-^\y - xf). 

Using (14. 4p and Lemma l3.H we have 



h <C2 \y - x\'^j{\y - x\)dy < C2 / 
J A Jbi 



\z\'^j{\z\)dz < 00. 



(4.4) 



(4.5) 



'B(0,R/4) 

For /3, we consider two cases separately: If < = 5^^{y) < fe(y), since v is decreasing, 

Hy)-K{y)<V{yd + R-'\y\^)-V{yd)= v{z)dz < R~'\y\Myd). (4.6) 

If yd = 5^^ (y) > doiy) and y £ E, using the fact that Soiy) is greater than or equal to the distance 
between y and the graph of ^pi and 

\y? 



yd-R+^/\W + {R-yd)^ 

we have 



I |2 I |2 

\y — x\ \y — x\ 



hx{y) - h{y) < 



Vd 



'R-V\y\^+iR-yd) 
Thus, by (|i:6]) - (f01) and Lemma [3Tl 



V|y|2 + (i?-y,)2 + (i?-y,) - 2{R-yd) " R 
v{z)dz < R~^\y - x|2 v{R - V|y|2 + {R- y^Y). 



(4.7) 



h < C2, I \y - x\'^j{\y - x\)dy < C3 \z\'^j{\z\)dz < 00. 

IE Jb{0,R/4) 
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We have proved 



\A{h-h^){x)\ <h+l2 + h<C4 



for some constant C4 = C4(i2, A) > 0. 

The estimate (j4.8p shows in particular that the limit 



lim / (Hy) - K{y))j{\y -x\)dy 

S-t-O J {y^D[jH+:\y~x\>£} 

exists and hence A{h — hx){x) is well defined. 

We now consider A(/i — hx). Note that for a.e. x G D n B{Q,R/4:), the second order partial 
derivatives of the function y — )• S^iy) exist at x. Without loss of generality we assume that x has 
been chosen so that the second order partial derivatives of the function y — )• S^iy) exist at x. Since 
hxiy) = V {(yd)'^) in CS, we have Ahx{x) = v'{xd). Moreover, since doiy) = yd for y = (0, Xd + e) 
when |e| is small, d"^ h{x) = v'{x^). Thus 



A(/.-M(x)=x: ''^S"'" u 



d-l 



i=l 



i=l 

dSpjy) 
dyi 



dyf 



E|-(^(^M.)) 



dSpjy) 
dyi 



\y=x 



\y=x 



+ v{6d{x)) 



dy. 



2 \y=^ 



(4.9) 



In the coordinate system CS, 



dyi 



\y=x 



and 



< — , i = 1, . . . , d — 1. 

- 3R' ' ' 



(4.10) 



Indeed, let e G M with |e| small, and x^^i := (0, . . . , e, . . . 0, Xd), i = 1, . . . , d — 1. Due to the uniform 
interior ball condition and the uniform exterior ball condition with radius R, we have 



R-^/^^T{R^^-Xd<5D{xe,i)-6D{x) < ^/e^ + {R + Xd? -R- 



Xd: 



SO 



-\SD{xe,i) - SDix)\ <- (^y^e^ + ^R-Xdf -{R- Xd)) V - (^V^^ + {R + Xdf -{R + Xd, 



V 



\^e^ + {R-Xd)^ + {R-Xd)J \^Je'^ + {R + XdY + {R + Xd)) 

which goes to zero as e — )• 0. The bound involving the second partial derivatives can be proved in a 
similar way using the elementary fact that ^ g^i^^ \ y=x = finie_s.o -^{^nixe^i) + ^D{x-e,i) — '^Soix))- 
Therefore, combining (j4.9p . (j4.10p and Lemma |3. 11 we have 



d-i 



\A{h-hx){x)\ <C5J2 



i=l 



dHD{y), 



dx 



2 \y=^ 



< 



C6- 



Using this, (jlT 
□ 



and linearity we get that A)h{x) is well defined and |(A+^)/i(x)| < cj. 
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We use C^iR"^) to denote the space of infinitely differentiable functions with compact support. 
Using the fact that A + ^ restricted to C^(M'^) coincides with the infinitesimal generator of the 
process X, we see that the following Dynkin's formula is true: for / G C^{W^) and any bounded 
open subset U of M"^, 

r (A + A)f{Xt)dt = E,.[/(X,J] - fix). (4.11) 
Jo 

Lemma 4.2 For every ri > and every a G (0, 1), there exists a positive constant c = c{ri,d,a) 
such that for any r G (0, ri] and any open sets U and D with B{0, ar) D D C U C D, we have 

F^{Xr^£D) < cr-'^K^[Tu], X £ DnB{0,ar/2). 

Proof. For fixed a G (0, 1), take a sequence of radial functions (j)m in C^{R'^) such that 1, 

0, \y\ < a/2 
i{y) = { 1, a < |y| < m + 1 
0, \y\ > m + 2, 



and that 11 V • • I -tv^ 



< Ci = Ci(a) < oo. Define (l)rn,r{y) = so that < (j)m,r < 1, 



4>m,r{y) = < 

Using (I4.12p . we see that 



0, |y| < ar/2 

1, ar < \y\ < r{m + 1) and sup ^ 
0, \y\ > r{m + 2), 



dyidyj 



< cir ^. 



(4.12) 



i,r{x + y) - 4'm,r{x) - (VfpmA^) ' y)'^B{o,i){y))J{y)dy 



< 



{\y\<^} 



.{x + y) -(t)m,r{^) - (V(/'m,r(a;) ' y)'^B{Q,l){y))J {y)dy 



+ 2 / J{:y)dy 

{\y\>i} 



<^ [ \y\^J{y)dy + 2 [ J{y)dy<^. 



(4.13) 



'{\y\<^} -^{ly^i} 
Now, by combining (I4.1ip . (14.120 and ()4.13p . we get that for any x £ Df] 5(0, ar/2), 

{Xru €{y€D:ar<\y\<{m+ l)r}) = E^, [0^,, (X,^) : Xr^ £ {y £ D : ar < \y\ < (m + l)r}] 



<E^ [^rr,,r (X,^)] =E^. 



{A + A)(t>ni,r{Xt)dt 



< ||(A + ^)(^„^,^||ooE^[7T/] < C4r-X-M. 



Therefore, since -8(0, ar) D D C U, 



(Xr^ eD)=¥^ (Xr^ e{yeD:ar< \y\}) 

= lim {Xr^ e{yeD:ar<\y\<{m+ l)r}) < C5 r-^E,^- 



□ 
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Define pq{x) := Xd — ^pq{x), where {x, Xd) are the coordinates of x in C5q. Note that for every 
Q e dD and x £ B{Q, R) Ci D we have 

(l + A^r^^PQix) < 5d{x) < Pq{x). (4.14) 

We define for ri , r2 > 

DQ{ri,r2) := {y G -D : n > pQ{y) > 0, \y\ < . 

Let Ri := R/{4:y/l + (1 + A)^). By Lemma EU V{5Dix) on the right-hand sides of (liTB ^ 
()4.16p can be replaced by 5d{x). The reason we prefer the forms below is that the function V will 
be used in the proof. 

Lemma 4.3 There are constants Aq > 2R^^ , kq G (0, 1) and c = c{R,A) > such that for every 
X> Xo, Q e dD and x G i:>Q(2~^(l + A)"^koA~\ kqA"^) with x = 0, 

P. e d) < cXV{6n{x)) (4.15) 

and 

< cX~^V{dD{x)). (4.16) 

Proof. Without loss of generality, we assume Q = 0. Let ip = ipQ : — M be the C^'^ function 
and CSf) be the coordinate system in the definition of C ' open set so that B{0,R)nD = {{y, ya) G 
B{0,R) in CSq : yd > ^(jj)} ■ Let p{y) := yd - ifiy) and D{a,b) := Do{a,b). 
Note that 

|yP = + IVdl"^ <r^ + ilVd -^iy)\ + \v{y)\? < (l + (l + for every y G D{r, r) . (4.17) 

By this and the definition of Ri, we have D{r, s) C D{Ri, Ri) C B{0, R/A) n D C B{0, R) n D for 
every r, s < i?i . 

Using Lemma l3.lt we can and will choose 6q G {0,Ri) small such that 

2r^ < T/((l + A2)-VV) for all r < ASq. 

Then, by ()4.14p . the subadditivity and monotonicity of V, for every A > 1 and every y G -6(0, R)riD 
with p{y) < 4A~^5o, we have 

2AV(y)' < V{X5Diy)) < (A + l)V{6Diy)) < 2XV{6D{y)). (4.18) 

Since A.ip{y) is well defined for a.e. y with respect to the (d— l)-dimensional Lebesgue measure, 
it follows that Ap{y) exists for a.e. y with respect to the d-dimensional Lebesgue measure. Using 
the fact that the derivative of a Lipschitz function is essentially bounded by its Lipschitz constant, 
we have for a.e. y G -6(0, R) Ci D that 

Ap{yf = A{yd - ^{y)f = 2(1 + \V^{y)\'') -2p{y)Aip{y) > 2(1 - p(y)||Ac^|U) • 

Choosing Jq ^ (0, -Ri) smaller if necessary we can get that 

Apiyf > 1 for a.e. y G -B(0, R) n D with p{y) < 25o. (4.19) 



Ex 



T-Dq{koA-i,A-1) 



18 



Let g{y) = g{y,yd) be a smooth function on M'^ with < g{y,yd) < 2, g{y,yd) < y% 



y\^\+y 



dyi 



< Cl, 



(4.20) 



and 



g{y) = < 



0, if - oo < yd < 0, or yd > 4 or \y\ > 2 

yj, if < yd < 1 and |y| < 1 

-(yd -2)2 + 2, if 1< yd < Sand |y| <1 

(yd -4)2, if 3 < yd < 4 and |y| < 1. 



Thus supp(y) C {|y| < 2, < yd < 4}. 

For A > 1, let gx{y) := gx{y,yd) ■■= y(A(5(7^y, A(5(7V(y)) so that 

supp(yA) C {|y| < 2A-^5o, < p{y) < 4A-^5o}. 

Then, since ■g{y)\ is essentiaUy bounded, using (|4.2U|) . we have 

d 



(4.21) 



E 



■9x{y)\ < C2>? a.e. y. 



(4.22) 



^9x{y) > A'^o"' for a.e. y G D{X"^6o, X~^So) 



dyidyj 

Note that, by the definition of g, gx{y) = \^5Q'^p{yY on L'(A"^(5o, A"Mo)- Thus, from (|4.19p we get 

(4.23) 

On the other hand, by (|4.22p we have 

/ {9x{y + z)- g\{y) - (Vgxiy) • z)1b{o,x-i){z))J{z) dz 

{gx{y + z) - gx{y) - i^gxiy) ■ z)iB{o,x-^){z))J{z)dz 



< 



+ 



{kl<A-i} 



{A-i<|z|<l} 



J{z)gx{y + z)dz+ ^ 



<C3A2 



{kl<A-i} 



|z|2j(z)dz + 2 



{l<kl} 



<kl<i} 
J(z) 



J(z)(iz yA(y) + 2 



{i<kl} 



J(z) 



+ 



{A-i<|z|<l} 



J{z)gx{y + z)dz + 



{A-i<|2|<l} 



J{z)dz gx{y). 
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Thus 



{gxiy + z)- g\{y) - {Vgx{y) ■ z)lB^o x-^){z))J{z)dz 

< C3 / \z\'^J{z)dz + 2X'~'^ / J{z)dz 

+A-2 / J{z)gx{y + z)dz + X-^l [ 

< C3 / \zfj{z)dz + 2X^'^ / J{z)dz 

Juz\<\-n Jn<\z\\ 



^<\z\<l} 



J{z)dz gxiy) 



+ 



{A-1<|2|<1} 



J{z)\z\'^gx{y + z)dz + 



{o<kl<i} 



\z\^J{z)dz gx{y). 



(4.24) 



We claim that for every A > 1 and y £ D{X ^Sq, X ^Sq), the function z — )■ gxiy + z) is supported 
inS(0,3A-i<5oV(4A)2 + l). 

Fix A > 1 and y G ^(A-iJo, A-^Jq) and suppose that z G S(0, 3A-^5o\/(4A)2 + 1)^ Then 
either li"! > 3A~"'^5o, or |z| < SA^-'^^o and \zd\ > 12X~^6oA. If |z| > 3A~"'^(5o, then clearly |y + z| > 
\z\ - \y\ > 3A-i(5o - X~^6o = 2X~^6o. Thus by (g^H), gx{y + z) = 0. Now assume \z\ < 3X-'^5o and 
l^dl > 12A"i(5oA. If < -12A-i(5oA, then ^A(y + z) = 0. If za > 12A-i5oA, we have 

p{y + z)>zd- \i>{y + z)\ > 12A-i5oA - k{\z\ + \y\) > X~^A{l26o - 3do - 6o) = 7X-^A6o. 

Thus by (I4.2ip . gx{y + z) = 0. The claim is proved. 

Using the above claim and the fact that gx{y) = ^^^'^piy)'^ on D{X^^5q,X~^5q), we have from 
(1424]) . that for y £ D{X~^6o, X'^Sq) 

X- 



{g\{y + z) - gx{y) - (Vgxiy) ■ z)lB(o,\-^)iz))Jiz)dz 

<C3 / \z\^J{z)dz + 2X~^ / J{z)dz 

J{\z\<x-^} Ai<kl} 



+ 



J{z)\z\''^dz + C4A^(5o ^yo(y) 



{A-i<|z|<lA3A-i5,)^(4A)2+l} 



<(C3 + 1) / 
ill; 



(lA|zp)J(z) d2 + 2A'" 



/ Jiz)dz + cA^6o^p{yf, (4.25) 
J{i<|^|} 



'{kl<3A-i5oV(4A)2+l} 
where C4 := 2"-*^ V /{o<|2|<i} {z]"^ J{z)dz. Define 

h{y) := V{6D{y))'i-Bio,R)nD(.y) and hx{y) := A/i(2/) - gx{y) 
Choose A* > 2 large such that for every A > A*, 



(C3 + 1) / 



{|^|<2A-i<5o^(4A)2+l} 



(lA|zp) J(z) dz+2X~'^ [ J{z) dz < 4-^(5(^2 
J{i<|^|} 



and T-^f^o — ^1' 



where Ci is the constant from Lemma |4.1[ Then by ()4.23p and (j4.25p . for every A > A* and a.e. 
y Gl?(A-i2-V^\,A-i(^o), 



(A + A)gx{y) > Agxiy) - \Agx{y)\ > A^^o" " ^-'>^'S^' - c^\%' p{yY > ^^6^' (4.26) 
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and 

(A + A)h,,{y) < A|(A + A)h{y)\ - (A + A)gx{y) < A(Ci - ^6^'') < -\x^6^^ . (4.27) 

Let 6^: := 2~^c^ ^^^^o a-nd / be a non-negative smooth radial function with compact support 
such that f{x) = for \x\ > 1 and f^a f{x)dx = 1. For k > I, define fk{x) = 2^'^f{2^x) and 

h^^\z):={fk*hx){z):= [ Uy)hx{z - y)dy . 

Let 

\ (k) (k) 

and consider large /c's such that S^'s are non-empty open sets. Since h)^ is in C^, Ah';;' is well 
defined everywhere. We claim that for every A > A* and k large enough, 

(A + A)h['^ < -l^'S^' on . (4.28) 

Indeed, for any x G and z G 5(0,2"'^'), when /c is large enough, it holds that x — z G 
D(A~^5*, A~^5o)- By the proof of Lemma l4.1l the following limit exists: 

lim / {hx{y - z) - hx{x - z)) i{\x - y\) dy 



£-S>0 



B{x,eY 



= lini / [hx{y') - hx{x - z)) j{\{x - z) - y'\) dy' = Ahx{x - z) . 

Moreover, by (j4.27p it holds that for every A > A,, {A + A)hx < -3A25-2 a.e. on D(A ^5*, A ^5o). 
Next, 

{h['\y)-h['\x))j{\x-y\)dy 



_ fk{z){hx{y - z)- hx{x - z)) dz ] j{\x - y\) dy 

\x-y\>£ 

fk{z) / {hxiy - z) - hx{x - z)) j{\x -y\)dy \ dz. 

5(0,2-*) \JB{x,eY J 

By letting e — t- and using the dominated convergence theorem, we get that for every A > A* and 
k large enough. 



iA + A)h['\x)= [ Mz){A + A)hx{x-z)dz<-]x\^ [ f,(z)dz = -]xH^^ 



By using Dynkin's formula (j4.1ip . the estimates (14.281) and the fact that h^-^^ are in C^(M.'^) 



X 

and by letting k ^ 00 we get for every A > A* and x G D{X^^5^, X'^^Sq) with x = 0, 

< -iA25o"'lE..[rB(A-i5„A-i5o)]- (4-29) 
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It is easy to see that hx > 0. In fact, if y G (-S(0, R) H Dy, then both h{y) and gx{y) are zero. If 
y e B{0,R)nD and p{y) > AX'^Sq, then gx{y) = 0. Finahy, if y G 5(0, i?) n and p{y) < AX'^Sq, 
then, since g{y) < y'^ by (|4.20p . we have from (j4.18p . 

hx{y) = XV{6D{y)) - g{X6^'y, XS^'p{y)) > XV{6D{y)) - X^p{y? > 0. 
Therefore, from (|4.29p . 

V{5d{x)) > ^A5o"'E-[^D{A-i5„A-i5o)]- (4-30) 

Since -B(0, (1 + A)"i(5^,A~i) nD C D{X~^6^, X-'^do), using LemmaSJand (|i30]) . we have that for 
every A > A* and x G B(0, 2~i(l + A)-M*A-^) with x = 0, 

P. G I?) < crX'E,,[ri,^x-^s„x-^So)] < cgAFl^^lx)). 

We have proved the lemma with Aq := X^^SZ^ . □ 



Lemma 4.4 There is a constant c = c{R,A) > such that for every X > Aq, k G (0, 1], Q G dD 
and X G DninX"^, X~^) with x = 0, 



^r,,^.-^ G Dq{2kX-\X-^) ) > cXV{5d{x)) 



(4.31) 



Proof. Fix A > Aq and k G (0, 1]. For simpHcity we denote Dq{kX ^, A ^) by I?. Further, let 

B = {ye D : pqiy) = kX~^ and \y\ < A~^} 

be the upper boundary of D. 

Let be the first time the Brownian motion W exits D and be the killed Brownian 

motion in D. Let Y = (Yj : t > 0) be the subordinate killed Brownian motion defined by Yj = W^^. 
Let C denote the lifetime of Y. Recall that u is the potential density of the subordinator S. It 
follows from |33l Corollary 4.4] that 



Thus, since u is deceasing, for any t > 0, 



¥,{Xrg £B) > E,. u{T^y,W^w G S,rf < tj > u{t)¥,{W^w ^ B,t^ < t) 



•d ' — — |_ -\ D 

u{t) W^{W^w G B) 



D 

> u{t) 



[W^w £ B) - t~^E,[T^ 



Now we use the following two estimates which are valid for the Brownian motion (for example, 
see |12t Lemma 3.4] with a = 0). There exist constants ci > and C2 > (independent of A > Aq) 
such that Fx(W^w G i?) > ciX6d{x) and E^; \t^] < C2X~^6d{x) ■ Then, by choosing to > so that 

D '- D ■' 

c\ - tQ^C2X~^ > ci - tQ^C2Xo > ci/2 =: C3, we get 

F^{X €B)> u{t){ci - C2t~^X~^)X6D{x) > C3u{to)X5D{x) . 



□ 
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5 Carleson estimate and Boundary Harnack principle 



In this section, we give the proof of the boundary Harnack principle for X. We first prove the 
Carleson estimate for X on Lipschitz open sets. 

We recall that an open set D in is said to be a Lipschitz open set if there exist a localization 
radius i22 > and a constant Ai > such that for every Q G dD, there exist a Lipschitz function 
tj) = tpQ : M'^~^ — )■ M satisfying V'(O) = 0, lipix) — ^(y)| < Ai|x — y|, and an orthonormal coordinate 
system CSq: y = (yi,... , yd-i,yd) =■ (y, Vd) with its origin at Q such that 

B{Q, R2)nD = {y= {y, y^) G 5(0, R2) in CSq : y^ > 

The pair (i?2, Ai) is called the characteristics of the Lipschitz open set D. Without loss of generality, 
we will assume throughout this section that i2i < 1. Note that a Lipschitz open set can be 
unbounded and disconnected. For Lipschitz open set D and every Q € dD and x G B{Q, R2) H 
we define 

Pq{x) := Xd - i'qix) , 

where {x,Xd) are the coordinates of x in CSq. 

The proof of the next lemma is similar to that of jl2|, Lemma 4.1]. 

Lemma 5.1 Let D d he a Lipschitz open set with characteristics {R2,Ai). There exists a 
constant 5 = 5(i?2, Ai) > such that for all Q G dD and x £ D with pq{x) < -R2/2, 

where t{x) := Tr)nB{x,2pQ{x)) = inf{i > : Xt i D r\ B{x, 2pq{x))}. 

Proof. Let := D n B{x,2pq{x)) and W^"" be the killed Brownian motion in Dj.. Here W 
denotes the Brownian motion in W^. As in the proof of Lemma [4. 4^ we define the subordinate killed 
Brownian motion Y = {Yt : t > 0) in by Yt := W^^{St). We will use C, to denote the lifetime 
of Y and let C^ := dD n B{x, 2pq{x)) and := mf{t > : Wt(^U}. 
Since, see M, G C^) > (Y^- G C^) = 



we have 



G D') > G C) > 



> u{t)F, W^w G C,,r]r < t > u{t) [FxiW^w G C, 



{T^>t)], t>0. (5.1) 



By the fact that D is a Lipschitz open set, there exists ci = ci{R2,Ai) > such that 

F^iW^w G a) > ci . 
(See the proof of |12i Lemma 4.1].) Since 



x(r^, >t)< 



t 



by using (j5.2|) and (jS.ip . we obtain that 



< 



2pqW)^ . (Pq(3^))^ , R2 

t '^^^ - '^T' 



(X,(,) G Z)^) > u{t) {^.{W^w^ G C,) - ¥^{t^^ > i)) > u{t) (^ci - C2^) 



(5.2) 



> ciu{to)/2 > 0, 
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where to = io{R2,-^i) > is chosen so that ci — C2Fl^/t > ci/2. The lemma is thus proved. □ 

Suppose that D is an open set and that U and V are bounded open sets with V C V C U and 
D n V ^ f/i. If / vanishes continuously on D"^ U, then by a finite covering argument, it is easy 
to see that / is bounded in an open neighborhood of dD n V. The proof of the next result is the 
same as that of [121 Lemma 4.2]. So we omit the proof. 

Lemma 5.2 Let D be an open set and U and V be bounded open sets with V C V C U and 
D r\V ^ %. Suppose f is a nonnegative function in M.'^ that is harmonic in D nU with respect to 
X and vanishes continuously on n U . Then f is regular harmonic in D OV with respect to X, 
i.e., 

/(x) =E,.[/(X,^^,,)] forallxeDnV. (5.3) 

Theorem 5.3 (Carleson estimate) Let D dW^ be a Lipschitz open set with the characteristics 
(i?2,Ai). Then there exists a positive constant A = A(i?2,Ai) such that for every Q G dD, 
< r < i?2/2, and any nonnegative function f in that is harmonic in Dn B[Q, r) with respect 
to X and vanishes continuously on n B{Q,r), we have 

f{x)<Af{xo) forxeDnB{Q,r/2), (5.4) 

where xq G DriB{Q,r) with pq{xq) = r/2. 

Proof. Since D is Lipschitz and r < ii2/2, by Proposition 12.21 and a standard chain argument, 
it suffices to prove ()5.4p for x G D n B{Q,r/12) and xq = Q. Without loss of generality, we may 
assume that f{xQ) = 1. In this proof, the constants 5,f3,r] and q's are always independent of r. 

Let v = z/(3)V2 where 1^(3) is the constant in (j2.8p with K = 3, choose < 7 < (z^~^ A (1 — 1^~^)) 
and let 

Bo{x) = DnB{x,2pQ{x)), Bi{x) = B{x,r^-^pQ{xy) 

and 

B2 = B{xo,pq{xo)/3) , B3 = B{xo,2pq{xo)/3). 
By Lemma l5.lt there exists 6 = 6{R2,Ai) > such that 

Px(X,so{.) ^D')>S, ^ e B{Q, r/4) . (5.5) 

By the Harnack inequality and a chain argument, there exists /3 > such that 

fix) <{pQ{x)/r)->'f{xo), xeDnB{Q,r/A). (5.6) 

In view of Lemma l5.2| / is regular harmonic in Bq{x) with respect to X. So 

fix) = E4fiXr,^J;X^^^^^^ G B^ix)] +E4fiX^^^J;X^^^^^^ ^ B,ix)] for x G B(Q,r/4). 

(5.7) 

We first show that there exists r/ > such that 

E4/(^-sow);^-so{.) ^ ^i(^)] ^ /(^o) ifx€DnBiQ,r/12) with pq(x) < r/r . (5.8) 
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Let rjo := 2 then, since ^ < I — u ^, for Pq{x) < rj^r, 

2pQ{x)<r^-^PQ{xy -2pq{x). 

Thus li X £ Dn B{Q,r /12) with pqix) < r]or, then \x - y\ < 2\z - y\ for z G Bo{x), y ^ Bi{x). 
Moreover, by the triangle inequahty, |x — y| < — z| + — y| < 1 + |z — y|. Thus we have by (j2.6p . 
(HZD, ([23]) and Lemma O 

=E.. / / j{\Xt-y\)f{y)dydt + E, / j{\Xt - y\)f{y) dy dt 

Jo J2>\y-x\>r^-ypQ{xp Jo J\y-x\>2 

<ciE,[TB,ix)]i [ j{\x-y\)f{y)dy+ f j{\x-y\)f{y)dy\ 

\J2>\y-x\>r'^-^PQ{x)-i J\y~'A>'i J 

<ciC2PQ{xf \ I j{\x - y\)f{y)dy 

\J\y~x\>r'^-ipQ{x)-i ,\y-xo\>2pQ(xQ)/'i 

+ [ j{\x-y\)f{y)dy] =: C3PQ{x)\h+l2). (5.9) 

J\y-'M<2pQixo)/3 J 

On the other hand, for z £ B2 and y ^ B3, we have \z — y\ < \z — xo\ + \xo — y\ < pQ{xo)/3+\xQ — y\ < 
2\xo — y\ and \z — y\ < \z — xq\ + \xo — y\ < 1 + \xo — y\- We have again by (|2.3p . (|2.6p . (|2.7p and 
Lemma 12.11 



f{xo) > f{Xr,J,Xr,.^^Bs 



>Exo/ "(/ j{\Xt-y\)f{y)dy+ f j{\Xt-y\)f{y)dy\dt 

; i|«-a;ol>2 y 



\J2>\y~xo\>2pQ{xo)/?. J\y-Ml 



\J2 



y\)f{y)dy + I j{\xQ-y\)f{y)dy\ 

2>\y-xo\>2pQ{xo)/'i J\y-'M>2 ) 

= C5Pq{xoY I j{\xo-y\)f{y)dy. (5.10) 

J\y-xo\>2pQ{xo)/3 

Suppose now that \y-x\> r^-^pQ{x)^ and x £ B{Q,r/A). Then 

\y - Xo\ < \y - x\ + r < \y - x\ + r'^pQ{xy^\y - x| < 2r'^ pQ{xy^\y - x\. 
Thus, using (|2.8p . we get for |x — y| < 2, 

j{\y - x\) < C7(/>Q(x)/r)~^^i(|y - xo\). (5.11) 
Now, using dlJ]), (EZl) and (l5TT]l . 



h < cr {PQ{x)/r)~^^3{\y-XQ\)f{y)dy 

J Ro/2>\y~x\>r^-i PQ{x)-i ,\y-xo\>2pQ{xo)/3 



+C8 / j(ko - y\)f{y)dy 

\y~-x\>Ro/2,\y~xo\>2pQ{xo)/'3 



< c<){{pQ{x)/r)-''' + i) / j(ko-2/l)/(y)dy 

J|j/-a;o|>2pg(xo)/3 

< C5-1c9/.q(xo)-2 ((pQ(x)/r)--^ + 1) /(xo) 

< 2c5-^C9(pQ(x)/r)-'^VQ(xo)-V(xo), (5.12) 
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where the second to last inequahty is due to (jS.lOp . 

If \y - xol < 2pq(xo)/3, then \y - x\ > \xo - Q\ - \x - Q\ - \y - xo\ > pq{xo)/6. This together 
with the Harnack inequahty imphes that 



h < cio 



j(k - y\)f{xo) dy < ciof{xo) / j{\x - y\) dy 

|t/-zo|<2pQ(a;o)/3 J \y-x\>pQ(x(,) /& 



Clofixo) / j{\z\)dz+ / j{\z\)dz 

\J Ro>\z\>pq{xo)/6 JRo<\z\ , 



< cio/(xo) / j{\z\)dz + cii . (5.13) 

\J Ro>\z\>pq{xo)/6 J 

Combining dO]), (fHT^ and (j^TIJ]) we obtain 

< ci2/(xo) (^pQ{x)'^{pQ{x)/ry^'''pQ{xo)~'^ 

{pQ{x)/rf{pQ{xo)/6f [ j{\z\) dz + {pQ{x)/TfA 

J B.o>\z\>pq(xo)/& ^ 

< ci3/(xo)((pQ(x)/r)2-7- + (pQ(x)/r)2( [ \z\''j{\z\)dz + l)\ 

\ ^JRo>\z\>pq{xo)/& ^ J 

< c,,f{x,){{pQ{x)/rf-^'' + {pQ{x)/rf) , (5.14) 
where we used the fact that pq{xq) = r/2. Since 2 — 71/ > 0, choose now rj G (0, ?/o) so that 

C14 {ri^-^" +rf) <l. 
Then for x G n B{Q, r/12) with Pq{x) < i]r, we have by (|5.14|) . 



We now prove the Carleson estimate (j5.4p for x G Dr\B[Q, r/12) by a method of contradiction. 
RecaU that /(xq) = 1. Suppose that there exists xi G Dr\B{Q, r /12) such that /(xi) > K > ?/~'^V 
(l+5~^), where K is a constant to be specified later. By (|5.6|) and the assumption /(xi) > K > r]^^, 
we have {pQ{xi)/r)^^ > /(xi) > K > r]~^, and hence Pq{xi) < rjr. By (|5.7p and (|5.8p . 



K < fix,) < E,., [fiXr^^^^^y,Xr^^^^^^ G i?i(xi) 

and hence 



+ 1, 



G i?i(xi)] > /(xi) - 1 > • 



In the last inequality of the display above we used the assumption that /(xi) > K > 1 + 6 ^. If 
K > 2/5/^, then L''^ n -Bi(xi) C -D^ n S(Q, r). By using the assumption that f = on DTi B{Q, r), 
we get from (j5.5p 



^xAfiXr,^,.,,),Xr,^^^^^^ G Bi(xi)] = E.J/(X.,^^^,^,),X,,^^^^, G i?i(xi)nZ)] 

< P,(X GD) sup /<(l-<5) sup /. 

' B^{x^) 
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Therefore, supB^(^^^-^ f > f{xi)/{{l + 6){1 — 5)), i.e., there exists a point X2 £ D such that 
\xi - X2\ < r^"'^ PQ{xiy and /(X2) > ^ fi^i) > . 

By induction, ii x^ £ D D B{Q,r /12) with /(xfc) > K/{1 - 5^)^-1 for k > 2, then there exists 
Xf^^i £ D such that 

\xk - Xk+i\ < r^~'^pQ{xkV and f{xk+i) > /(^^O > ^-^^ -^^2)fc ^ " ^^'^^^ 

From dHj]) and (l5J5]l it follows that PQ{xk)/r < (1 - ,^2){fc-i)//3j;^-i//3^ foj. g^ery /c > 1. Therefore, 

k—l 00 

< \xi- Q\+^\xj+i- Xj\ < — + ^r^"'^pQ(xj)^ 

00 00 

<— + r^"^ ^(1 - 52)0-l)7//3^-7//3^7 = ^ + ^l-7^7_^-7//3 _ j2yV/9 

■"■^ i=i """^ 3=0 



12 l_(l_52)7//3- 

Choose 

K = rjy{l + 5-^) V 12'^/^(1 - (1 - 52)7//3)-/3/7. 

Then K-t//^ (1 _ (1 _ < 1/12, and hence Xk £ D r\ B{Q,r/6) for every A; > 1. Since 

limk-^oo f{xk) = +00, this contradicts the fact that / is bounded on B{Q,r /2). This contradiction 
shows that f{x) < K for every x £ D r\ B{Q, r /12). This completes the proof of the theorem. □ 



Proof of Theorem [m . We recall that Ri = R/{4^yi + (1 + A)^) and Aq > 2i?j"^ and kq £ (0, 1) 
are the constants in the statement of Lemma 14.31 

Since D is a C^'^ open set and r < i?, by the Harnack inequality and a standard chain argument, 
it suffices to prove (jl.4|) for x,y £ DD B{Q, 2~^rKoXQ^). In this proof, the constants r] and Cj's are 
always independent of r. 

For any r £ {0,R] and x £ D D B{Q,2~^rKoXQ^), let Qx be the point Qx £ dD so that 
l^; — Qx\ = Sd{x) and let xq := Qx + |(x — Qx)/\x — Qx\- We choose a C^'^-function 99 : M'^"-'^ — )• M 
satisfying (^(0) = 0, Vip{0) = (0, ...,0), \\Vip\\oo < A, \Vip{y) - Vip{z)\ < A\y - z\, and an 
orthonormal coordinate system CS with its origin at Qx such that 

B{Qx, R)nD = {y= {y, y^) £ B{0, R) in CS:yd> ^{y)}- 

In the coordinate system CS we have x = and xq = (0, r/8). For any bi, 62 > 0, we define 

D{bi,b2) := {y = {y,yd) inCS :0 <yd- (p{y) < biVKoXo^, \y\ < b2rXQ^} . 

It is easy to see that D{2, 2) C DnB{Q,r/2). In fact, since A > 1 and R<1, for every z £ D{2,2), 

\z - Q\ < \Q - x\ + \x - Qx\ + \Qx - z\ < \Q - x\ + \x - Qx\ + \zd - (piz)\ + 1^(2)1 
< rAo ^((1 + A) + 4) < 2-Vi?((l + A) + 4)/(4Vl + (1 + A)2) < ^. 
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Thus if / is a nonnegative function on M'^ that is harmonic in D n B{Q,r) with respect to X and 
vanishes continuously in D^r]B{Q, r), then, by Lemma [5.2| / is regular harmonic in DriB{Q,r/2) 
with respect to X, hence also in D(2, 2). Thus by the Harnack inequality, we have 



f {-^Tori.!)) 



G -0(2,1) 



(5.16) 



> ci/(xo)P4X,^(,,,) G D{2, 1)1 > C2f{xo)6D{x)/r. 



In the last inequality above we have used (I4.3ip . 

Let w = (0, tAq ^Ko4). Then it is easy to see that there exists a constant rj = 7?(A, 6q) € (0, 1/4) 
such that B{w,r]rXo^Ko) G D{1, 1). By ([M]), dlZD, and Lemma[2Tl 



/(u;) > E^ /(X,^(,^^j);X,^(^^^j ^Z)(2,2 
> C3E„, [r 
Hence by HMD, dZZl), (HIS]), 



E„ 



Rrf\D(2,2) 



/(y)i(l^t - yD^iydi 



BKr^rA-^Ko)] / f iv) J il'^ - vl) dy > c^r^ 



R<'\Z)(2,2) 



f{y)j{\w - y\)dy. 



/ ^I)(2,2) 



E. 



Td{1,1) 



< C5E^[r£,(i^i)] / 

< C6fe(a;)r 



\D(2,2) 



f{ymXt-y\)dydt 



M.<i\D{2,2) 



f{y)i{\w - y\)dy 



f{y)j{\w-y\)dy < 2^M^f{w). 

I8'*\D(2,2) C4 r 

On the other hand, by the Harnack inequality and the Carleson estimate, we have 



E, 



/ ; G -0(2,2) < C7/(xo)P. (Xr^,,,,^ G D{2,2)\ < cs f{xo)SD{x)/r. 



In the last inequality above we have used (|4.15p . Combining the two inequalities above, we get 



/(^) 



/ ei)(2,2) 



+ E, 



/ P^^z,a,i) ;^-z.<M, ^^(2>2) (5.17) 



< %^{x)fixo) + ^^^^f{w) < ^5D{x){f{xo) + f{w)) < ^6D{x)fixo). 
r c^r r r 

In the last inequality above we have used the Harnack inequality. 

From (j5.16p - (|5.17p . we have that for every x,y E D Ci B{Q, 2~^rKQXQ^), 

f{x) ^ cio Spjx) 
f{y) ~ C2 Sniy)' 



which proves the theorem. 



□ 
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6 Counterexample 



In this section, we present an example of a (bounded) C^'^ domain (open and connected) D on 
which the boundary Harnack principle for the independent sum of a Brownian motion and a finite 
range rotationally invariant Levy process fails, even for regular harmonic function vanishing on D^. 
A similar example appears in \20\ Section 6] for the case of truncated stable process. 

Suppose that Z is a rotationally invariant Levy process whose Levy measure has a density 
J{x) = j{\x\) with j(r) = for all r > 1. Suppose that Z is independent of the Brownian motion 
W. We will consider the process Y = W + Z . For any Borel sets U and V in R with V C U , we 
have 



■,{Y^y€V)=E, ji\Yt-z\)^Y^^,\^i}i\Yt-z\)dzdt xGU, (6.1) 



V 



where := inf{t >0:Yt^U}. 

We consider the bounded domain in 



D := (-100,100)" \ (^(-100,49]'^"^ x [-1/2, 0]j . 

Suppose that the (not necessarily scale invariant) boundary Harnack principle is true for Y on 
D at the origin for regular harmonic function vanishing on D^, i.e., there exist constants Ri > 
and Ml > 1 such that for any r < i?i and any nonnegative functions u, v on M.'^ which are regular 
harmonic with respect to y in D n B{0, Mir) and vanish in D^, we have 

-rT<Cr-r^ for any G L>nS(0,r), (6.2) 

v{x) v{y) 

where Cr = Cr{D) > is independent of the harmonic functions u and v. Choose an ri < Ri with 
Miri < 1/2 and let A := (0, ^^i). We define a function v by 



vix) ■.= F^[Y^Y £{yeD;yd>0} 

By definition v is regular harmonic in Dr]B{0, Miri) with respect to Y and vanishes in D^. Applying 
the function v above to (j6.2p . we get a Carleson type estimate at 0, i.e., for any nonnegative function 
u which is regular harmonic with respect to y in Z) H B(0, Miri) and vanishes in we have 

v(A) 

u(A) > c-^--^u(x) > c-^v(A)u(x) = ciu(x), x G £> n 5(0, n), (6.3) 
^ v{x) ^ 

where ci = c~^v{A) > 0. We will construct a bounded positive function u on M'^ which is regular 
harmonic with respect to Y m D f] 5(0, Mir) and vanishes in for which (|6.3p fails. 
For n > 1, we put 

Cn ■■= {ix,Xd)£D: |x|<2-"-3ri, X2< -1 + 2^^2}, 

Dn ■■= {{y, Vd) & D : yd>0, |x - y| < 1 for some x G Cn} ■ 

It is easy to see that 

a: C {(y, yd) : \y\ < (2-"-^ + 2-("-i)/2)ri, < yd < 2-'^r'\) C 5(0, ri) n L», for n > 2. (6.4) 
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Indeed, for any y G Dn, we have yd S [0,2 "rf] and |y — j;| < 1 for some x S C„. If \y\ > 
^2-n-3 _^ 2-("-i)/2)^^^ yd > and X G C„, then 

\x - > + (|y| - \x\f > (1 - 2-"r?)2 + 2-("-^)r? > 1. 

Thus, in this case y ^ 

For any n, let be the first hitting time of Dn by the process Y. By ([6 

-I' ^ T^i" \ . TO . f^Y 



y'Dr\B{o,Miri) > j (^''"z)nB(o,Afiri) > ^{0} j ~ ^' as n -)> oo. 

Fix no > 2 large so that 



and define 



ff^A(T^nB(o,A/in)>^^no) < I (6-5) 
n(x) := (y^y e C„„ 

\ DnB{0,Miri) 

u is a nonnegative bounded function which is regular harmonic in D H B{0, A/iri) with respect to 
Y and vanishes in D^. It also vanishes continuously on dD n -6(0, Miri). Note that by (|6.ip . 



^CnB(0,A/iri) ° -Uni((,U,A7irij — Un,^ J ^ y ^CnS(0, A/j ri ) \D„g ° > 

Thus by the strong Markov property, 



, 

Y ^ rrY 



^ °^ \xeD„o ) 2 yxgDnB(0,ri) y 

In the last inequality above, we have used (|6.4p - (|6.5p . But by ()6.3p . > ci sup^g^ip^^o.n) ""C^)' 

which gives a contradiction. Thus the boundary Harnack principle is not true for D at the origin. 

By smoothing off the corners of D, we can easily construct a bounded C^'^ domain on which 
the boundary Harnack principle for Y fails at 0. 



7 Proofs of Theorems 11.41 and 11.51 

As already said in the introduction, once the boundary Harnack principle has been established, 
the proofs of Theorems 11.41 and 11.51 are similar to the corresponding proofs in [13] for the operator 
A + a" A"/^. In fact, the proof are even simpler, because [13] strives for uniformity in the weight a. 

The proof of Theorem ll.4l in the case (i > 3 is by now quite standard. Once the interior estimates 
are established, the full estimates in connected C^'^ open sets follow from the boundary Harnack 
principle by the method developed by Bogdan [3] and Hansen [19]. For the operator A + aA"^/^ 
this is accomplished in [131 Section 3] . In the present setting the proof from |13] carries over almost 
verbatim. In several places in |13j one refers to the form of the Levy density, but in fact, the form 
of the Levy density is only used to establish uniformity in the weight a. 
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When d = 2, the above method ceases to work due to the nature of the logarithmic potential 
associated with the Laplacian. The proof in [131 Section 4] for the operator A + aA"/^ uses a ca- 
pacitary argument to derive the interior upper bound estimate for the Green function. By a scaling 
consideration and applying the boundary Harnack principle, one gets sharp Green function upper 
bound estimates. For the lower bound estimates, [13] compares the process with the subordinate 
killed Brownian motion when D is connected, and then extend it to general bounded C^'^ by using 
the jumping structure of the process. In the present setting, the proof of the lower bound is exactly 
the same as in [13] (see proofs of Theorems 4.2 and 4.4). The proof of the upper bound is essentially 
the same as the one in [T3], except that one has to make several minor modifications. Lemma 4.5 
in [13] should be replaced by the following statement: There exists c > such that for any L > 0, 



This is proved in the same way as [13^ Lemma 4.5] by using the explicit formula for the Green 
function of the ball B{0, L) C M^: 



The statement of Lemma 4.6 in [13] should be changed to: There exists c > such that for any 
L > and bounded open set D in containing -6(0, L) and any x G -6(0, ^) 



c 



"""^"'"^ - CapMi?(0,N/2)) '^-^"^™^"" 

(we refer to [13] for all unexplained notation). Next, Corollary 4.7 in [13] should be replaced by 
the statement: There exists c > such that for any L > and any x G -8(0, 3-L/4) 

Gb(o,l){x,0) < clog(L/|x|). 

Finally, the last change is in the proof of Lemma 4.8 in ^3] which uses a scaling argument. This 
in our setting can be circumvented by using the modified statement of [1'6\ Lemma 4.6]. The rest 
of the proof remains exactly the same. 

The proof of Theorem 11.51 is also quite standard (see [2} 113 1 I^T ] I22j). In the current setting we 
follow step-by-step the proof of the corresponding result in |13l Section 6]. The main difference 
is that [13] uses the explicit form of the Levy density for the operator A + aA"/^ which is 
c{a, d, a)r~'^~°' . This Levy density is now replaced by j, and it suffice to use properties (j2.6p and 
()2.7p to carry over all arguments. The reader can also compare with [221 Section 6] where the 
Martin boundary was identified with the Euclidean boundary for purely discontinuous processes 
whose jumping kernel satisfies (|2.6p and (|2.7p . 
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